REGULARITY CONDITIONS IN THE CLT FOR LINEAR 
EIGENVALUE STATISTICS OF WIGNER MATRICES 
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Abstract. We show that the variance of centred linear statistics of eigenvalues of GUE 
matrices remains bounded for large n for some classes of test functions less regular than 
Lipschitz functions. This observation is suggested by the limiting form of the variance 
(which has previously been computed explicitly), but it does not seem to appear in the 
literature. We combine this fact with comparison techniques following Tao-Vu and Erdos, 
Yau, et al. and a Littlewood-Paley type decomposition to extend the central limit theorem 
for linear eigenvalue statistics to functions in the Holder class (7 1 / 2 + e j n the case of 
matrices of Gaussian convolution type. We also give a variance bound which implies the 
CLT for test functions in the Sobolev space H +e and C 1_e for general Wigner matrices 
satisfying moment conditions. Previous results on the CLT impose the existence and 
continuity of at least one classical derivative. 



1. Introduction 
We are concerned with the fluctuations of linear eigenvalue statistics 

n 

i=i 

where Xj, 1 < j < n are the eigenvalues of a Hermitian Wigner matrix. The centred 
random variable 

(1) A/>]:=AU^]-E(A/-„M) 

exhibits strong cancellation, and it is known that for sufficiently smooth test functions, it 
converges to a normal random variable. This result is originally due to V. Girko, see [19J. 
In contrast, sums of n centred i.i.d. random variables require a normalization of order n~2 
to obtain a non-trivial Gaussian limit. 

Some degree of smoothness is required of ip for the asymptotic normality to hold without 
normalization, and the variance to remain bounded. Indeed, J. Gustafsson [T7] has shown 
that for GUE and for y in the interior of the support of the limiting spectral density, and 

ifyix) = l[y i0o) (x), 
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i.e. J\f n [(p y ] is the number of eigenvalues greater than y, the variance has logarithmic 
asymptotic behaviour: 



\ar(NnW\) = I ^2 ] lo e"- 

Moreover, Af°[ip y ]/ ^/~Var(M n [ l Py)] converges to a standard normal random variable. S. 
Dallaporta and V. Vu [9] have recently extended this result to Wigner matrices whose 
entries have sub-exponential tails. 

It is natural to ask what is the minimal degree of smoothness required of (p for the 
variance to remain bounded, and for N° [ip] to converge in distribution. The present work 
is motivated by this question. 

In recent years, several authors have derived central limit theorems for centred linear 
eigenvalue statistics of Wigner matrices with various hypotheses, see [I], [5], [7], [6], [9], 
[20] . [21], [26], [29], [31], [32]. In all these works, when general test functions are considered, 
these are required to be at least C 1 (IR) (differentiable with a continuous derivative) and 
often much smoother. For example, A. Lytova and L. Pastur [26j have proved that for GOE 
matrices and cp 6 C 1 (M) satisfying suitable decay conditions at infinity, N°[tp\ converges 
to a Gaussian random variable with mean zero and variance 



The variance Var(J\f n [p]) is controlled by a Poincare-type bound specific to the Gaussian 
case. Under a Lindeberg-type condition on the fourth moments, Lytova and Pastur extend 
their result (by comparison to the Gaussian case) to general real Wigner matrices assuming 
(p £ C 5 (M). For general symmetric Wigner matrices, there is an additional term in the 
expression for limiting variance on the right side of Q, which vanishes if 

E|u>jj| 4 = 3B\wij\ 2 , 

where Wij = Wij are the entries of the Wigner matrix. One expects the fluctuations of 
linear statistics to be sensitive to the third and fourth moments of the entry distribution. 

G. Anderson and O. Zeitouni [lj have developed a central limit theorem for linear sta- 
tistics of C 1 (M) test functions for a large class of random matrices with independent, but 
not necessarily identically distributed entries which includes Wigner and Wishart matrices 
with entries satisfying a Poincare inequality. 

Sobolev spaces offer an alternative to classical derivatives as a way to measure smooth- 
ness. For s > 0, the (inhomogeneous) L 2 -Sobolev space H S (K) is defined as the closure of 
the space of Schwartz functions S in the norm: 



2 / f-t i ltl\2s 



i + ki) 2s i/(oi 2 de, 
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where / is the Fourier transform of /. K. Johansson |21j considered Hermitian matrix 
models defined by their density with respect to Lebesgue measure on the entries: 

d/i„(x) = ^exp(-/?iVTW(x))dx, 
dx = dtftxijdQxij dx„. 

i<j i 

It is well-known that for such matrices, the joint density on ~K N of the eigenvalues has the 
form 



— ex P ( -N £ V(\j) + log |Ai - Xj\ J 



and that the classical GUE ensemble corresponds to /3 = 2, V = 2x 2 . Johansson obtains 
the central limit theorem for linear statistics for test functions of class H 2+e (IR) when (3 = 2 
and H 17 / 2+e (R) when (3 ^ 2^ Based on an idea in [21], M. Shcherbina [29| (see also [30]) 
proved a variance bound which allowed her to extend the central limit theorem for linear 
statistics of real Wigner matrices to functions in the Sobolev space ff 3 / 2+e (R), for e > 
and under weak assumptions on the moments of the entries. Note that such functions 
are in C 1 (M) by the Sobolev embedding. The key estimate in Shcherbina's work is the 
following bound for the variance of Af n [<p] (see [29], Prop. 1): 

(3) Var(AT n M) < C s \\cp\\ 2 Hsm / e'V^ 1 / Var(TrG(£ + irj)) dEdrj, 

JO J-oc 



where 



G(E + i?]) = G(z) ' 



H - z 

is the resolvent matrix. It is the formula on the right side in ^ which motivated the 
Littlewood-Paley approach explained in Section [3| Equation ^ shows that it suffices to 
obtain a bound for the variance of TrG(z) of order 7]~ 2s to conclude that the variance of 
linear statistics can be bounded in terms of the H s+e norm of the test function for e > 0. 
(The integration in x turns out to be harmless). 

In their work on universality for Wigner matrices, L. Erdos, H. T. Yau and their col- 
laborators have obtained large deviation bounds of order r]~ 2 n 2 logn for the trace of the 
resolvent TrG(z) for r/ S> n _1 (see the survey paper [llj). We use their methods to show 
that 

(4) Var(TrG(z)) < C e rT 2 ~ e , < 5. 

Our result for linear statistics of functions in H 1+e follows from this estimate: 



In [21], Remark 2.5, K. Johansson explains that the restrictions if € H 2+e and <p G H 17 ^ 2+c are technical 
and that the "correct condition should be the finiteness" of the limiting expression for the variance. 
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Theorem 1. Let tp £ H 1+e (R), where e > is arbitrary. Define the random variable 

n 

where Xj , 1 < j < n are the eigenvalues of a Hermitian Wigner matrix whose entries satisfy 
the following condition ("Condition CO"): 

(5) P(Kj| > t°) < e~\ c>0. 
There is a constant C > such that 

(6) Var(W n [c/?]) < C\\ip\\ 2 Hl+HR) . 

In particular, M^[(p\ converges in distribution to a random variable with variance given by 

f2 r 2 



2 J -2 



4^ 2 .Lo V a: - 2/ / ^/4^2 ^73^2 



where 



k 4 f f 2 , , 2-x 2 , \ 2 W2-2{[ 2 . . 2-x 2 ^ 2 

K4 = E|?«jj| 4 — 3 
w 2 = E|u;jj| 2 . 

For GUE matrices, the variance of the limiting random variable becomes (see [5], or [28] 
Ch. 5): 



This integral is finite for a- Holder functions <£> of order a > 1/2 (see Section [5]). For 
< a < 1, the space C a (IR) of a- Holder functions consists of (continuous) functions such 
that the norm: 

IMU = IMIl°°(r) + Mc*«(r) 

Mc^sup '^-y 

x ¥ = y I*" 2/1 

is finite. Since we are only concerned with functions defined on the real line, we use the 
notation C a = C a (R) throughout. 

One expects that also for finite n the variance of Var(A/" n [<p]) can be bounded in terms 
of the C 1//2+<E norm of 92. This forms the basis of our second result: 

Theorem 2. Let <p £ C l l 2+e for e > arbitrary. Define the random variable 

n 

(8) Mn[<p]=Y,<p(\j), 
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where Xj are the eigenvalues of a GUE matrix of size n. There exists a function C{K) > ; 
bounded on compact sets, such that for all ip € C l l 2+t , \\(p\\ c i/ 2 +e < K 

(9) Var(AT n M) <C{K) |M|*, 1/a+e . 

Moreover, the centred random variable 

CM = MiM - E(jV" n [^]) 

converges in distribution to a normal random variable of mean zero and variance given by 
&■ 

By a moment matching argument and a saddle point analysis, we obtain a central limit 
theorem for more general Hermitian matrices, of the type considered by K. Johansson in 
[22]. 

Theorem 3. Let W = {wij)i<i,j<n be a Hermitian Wigner matrix whose entries satisfy: 

P(Kj| > t c ) < e~\ c> 0, 

and such that the first five moments of Wij match those of the GUE. Form the Wigner 
matrix M 
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+ V)i an d define the random variable Af n [ip] as in ^8), 



but with Xj 



being the eigenvalues of M . Then M^if] converges in distribution to a mean zero normal 
random variable with variance given by 

1 (<p{x)-<p{y)\ 2 2-xy 



dxdy. 



The variance normalization in the previous theorem is of no particular significance. We 
have chosen it as a matter of convenience, for the calculations appearing in Sections [6] and 

m 

In the final section, we derive the following central limit theorem for linear eigenvalue 
statistics for general Wigner matrices under matching moment conditions for regularity 
below one bounded derivative: 

Theorem 4. Let 

H = -^=W 
in 



be a Wigner matrix such that the entries of W satisfy condition (CO) and such that the 
distribution of its entries matches that of GUE up to five moments. Then for e > 
sufficiently small, and any function ip £ C 1- % the centred random variable 

■CM = N n [p] - E(A/"nM) 
converges in distribution to a normal random variable of mean zero and variance given by 
&■ 

The results above do not provide sharp characterizations of the functions for which the 
variance remains bounded and the CLT holds. Indeed, we expect the variance to remain 
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asymptotically bounded for all functions <p such that ([7]) is finite, under weaker assumptions 
on the matrices than the ones we make. The integral 

(x - y) z 

is equal, up to a constant factor, to the homogeneous i7 1 / 2 (IR) norm: 

n/y = / ien/(e)i 2 de 

For functions supported away from the spectral edges ±1, the regularity threshold for 
asymptotic normality to hold without normalization is likely to be H 1 / 2 . This is consistent 
with alternate expressions for the variance found in [5] and [21], as well as the logarithmic 
divergence of the variance in case <p = lj, the indicator function of an interval / in the 
interior of the limiting spectrum. Such functions barely fail to be in H 1 / 2 . Note in that in 
the case of Haar-distributed matrices on the unitary group, P. Diaconis and S. Evans [10] 
have obtained a sharp (with respect to the regularity of the test functions) central limit 
theorem. They find that the variance of the analogue of Af n [<p] grows like the nth partial 
sum of the appropriate H 1 / 2 norm. However, their techniques rely on exact computations 
of matrix traces that are not available in the case of Wigner matrices. 



Acknowledgements. The first author would like to thank David Renfrew for pointing 
out M. Shcherbina's result [25], and for discussions that provided the initial motivation for 
the present work during the 2011 Prague Summer School on Mathematical Physics. We 
would also like to thank Laszlo Erdos for helpful comments on the variance estimates ([3]) 
and (fill. 



2. Approximation by Poisson Integrals 

To derive the central limit theorem for functions with limited smoothness, we bound the 
variance of A/^[y] by the norm of ip in a function space describing the regularity of <p. The 
central limit theorem then follows by a simple approximation argument. In what follows, 
we use the next lemma, adapted from [25], Proposition 3.2.5 (see also [26J, p. 1793). 

Lemma 2.1. Suppose X is a normed vector subspace of Cfe(lR), i.e. continuous functions 
bounded on M, such that for \\<p\\x < K, we have the bound 

V n (<p) := Var(.A/*„M) < C(K)\\<pf x , 

where C : M + — > M + is bounded on compact subsets. If there exists a dense subspace Y C X 
such that for ip £Y , we have 

Z n [<p] := E(e^M) = exp(-xV(^)/2) + o(l), 

for a continuous quadratic functional V on X , then the central limit theorem is valid for 
all ip £ X . 
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Proof. Let {ipk} C Y converge to cp E X. In particular, HV'fcllx is a bounded sequence. 
Then: 

\Z n [^ k }(x) - Z n [<p](x)\ < \x\B [\N°[i> k ]-KW\ 

< \x\Vsx(N n [<p - ^ k )) 1/2 

< \x\\\ip - ipk\\x- 

The result follows from the continuity of V and the Levy continuity theorem. □ 

In her proof of the CLT for linear statistics for H 3 / 2+e functions |29j . M. Shcherbina 
used approximation by functions of the form: 

where P^ is the Poisson kernel 

1 7] 

p y( x ) = — ^ — 2' V > °- 
tt rj z + x z 

When / € L P (IR) for l<p<oo, P r? */is real analytic, so the validity of the CLT for 
such functions follows from work anterior to [29] . see [31], However, since P^{x) = 3 s x ^ , 
the CLT for Poisson integrals has a direct proof based on a calculation involving resolvents 
(see the proof of Theorem 1 in |29j). 

We use the following standard approximation results for Poisson integrals: 

Lemma 2.2. Let f £ H S (M), for < s < oo. Then the Poisson integral P v * f converges 
in the H s norm to f for r\ — > 0. 

Lemma 2.3. Let f £ C a (M) for some < a < 1, and f3 > be less than a. The Poisson 
integral P-q* f converges to f in the norm for 77 — >• 0. 

The first result is a simple consequence of the explicit form of the Fourier transform of 
P v *f- 

From this, we have: 

\\Pn*f- f\\h = /(I + I£D 2S (1 " e- m ) 2 \f(0\ 2 d£, 
and the Lemma 2.2 follows from the dominated convergence theorem. 



We indicate how to derive Lemma |2.3| We begin with the L°° part of the norm: 
|(P„ */)(*) -/(*)!< / \P v (t)\.\f(x-t)-f(x)\dt 



< [ \P v (t)\\t\ a dt. 

JR 



< n a 

r^j '/ 1 
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so P-q * / converges uniformly to / as r\ tends to 0. The convergence of the seminorm 
[P y * f] a can be proved using Littlewood-Paley techniques, or using the characterization 
found in [34, Section V.4]; for / G C a , the seminorm [f] a is equivalent to the least M such 
that 

||<VWII°o<M7T 1+a 

Thus 

[P v *f-fh<v 1 ^\\d v P v *f\\oc 

The convergence follows from this. 

3. Littlewood-Paley Theory 

In this section, we derive a representation for the variance of A/" n [<£>] based on a decom- 
position of ip of Littlewood-Paley type. 

Theorem 5. Let < s < oo and < p < a < 1. Then if tp G H s or ip G C a , ip has a 
representation: 

oo 

= ( P 2-fc *9k)(x), 

k=-l 

where the sum on the right is convergent in H s if ip G H s , and in if ip G C a . T/ie 
functions gk, — 1 < A; < oo are smooth and such that, for some constants C s and C a ^: 

oo 

E 2fcsii i|2 ^ || i|2 
2 \\9k\\ L 2 < C s \\ip\\ HS 

k=-l 

and 

oo 
fc=-l 

Corollary 6. Lei (/? G 5(IR), i/ie cZass of Schwartz functions, i.e. C°° functions on M such 
that sup xg]R \ x m f( k \x)\ < C m k < oo for m, k < oo. The variance of M n [ip] can be written 
as 

oo r r 

(10) Var(AA n M)= ^ / g k (t)gJs)Cov(M n [P 2 - k (- - t)],M n [P 2 -i(- - s)]) dtds. 
k,i=-i J R J R 



The sum on the right side of (10) is absolutely convergent for ip G S. 

As announced above, the main tool in the proof of Theorem [5] will be Littlewood-Paley 
theory. A good general reference on this is [2j Ch. 2]. The dyadic frequency decomposition 
was introduced in the context of Fourier series by J.E. Littlewood and R.E.A.C. Paley [23], 
|24j . [25]. E. M. Stein developed extensions of the theory to functions on higher dimensional 
spaces M d using various Littlewood-Paley functions based on the Poisson kernel (see 



Chapter IV]). Here we will use the smooth dyadic frequency decomposition presented in 
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[2] (see also [8]), and the description of the functions spaces C a and H s in terms of this 
decomposition. We denote by 

7(6 = / e-^f(x)dx 

the Fourier transform of / € S', the space of tempered distributions. 

Given a function ip £ S, its inhomogeneous Littlewood-Paley decomposition is the sum 

k>-l 

with 

<£>_i = h* {p 
ip k = 2 k u(2 k -)*ip, k>0, 

where < h < 1 and < Q < 1 are even functions with compact support in the ball 
< 3/4} and the "annulus" C = {3/4 < |f| < 8/3}, respectively, such that: 

k>0 

o< no 2 + Y J Q2 ( 2 ~ k o ^ L 

fc>0 

for Note that iff. is frequency- localized in an annulus of size 2 fc , i.e. 

supp <p~k C 2 k C. 

The first component ip_i has its Fourier support contained in a ball 

supp ipZi C M < 3/4}. 

For the existence of this decomposition and the corresponding dyadic partition of frequency 
space, see [2], Proposition 2.10. 

The following theorem characterizes H s (R) and C a (R) in terms of the Littlewood-Paley 
decomposition (see [2], Section 2.7): 

Theorem 7. T/ie inhomogeneous Sobolev space H S (R) consists of those tempered distri- 
butions such that 

II „l|2 2fcsii „ 1 1 2 , _ 

IMbf 2 ■= 2^ 2 H^felli 2 < °°- 
k 

Moreover, the norms \\ ■ ||b| 2 and \\ ■ \\h s are equivalent. The space C a (M) of a-Holder 
continuous functions is the space of distributions such 

\\<p\\b°> x ■= sup2 afc ||93 fc || L oo < oo. 

k 

The norms \\ ■ and \\ ■ ||c a are equivalent. 
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We turn to the proof of Theorem [5] and its corollary. Given the support properties of 
tpk, we can define gk by 

We have: 

\\9k\\iA> ^ IIVfellLP, P = 2,oo. 
Indeed, for k > 0, let i[)(2~ k •) be a function in C 00 which is 1 on the support of the 
Littlewood-Paley cutoff uj(2~ k •). Then: 

J 7-1 denotes the inverse Fourier transform. By Young's inequality ([2], Lemma 1.4) 

||9 fc ||Lp<||^ 1 (e 2 - fe l^(2-^))llLi||^ll^, 
it suffices to obtain an L 1 bound for 

F-\e 2 - k \^{2- k i)) = /e**e 2 ~ fc l^(2- fc 0de. 



Let us compute: 

1 1 J e^e 2 ' k ^(2- k 0d^ 
Changing variables s = 2~ k £, we have: 

e ix2ks e^iP(s)2 k ds 



dx. 



dx. 



Another change of variables (in the outer integral) t = 2 k x gives: 



? its e^iP(s)ds 



dt. 



The integrand is 0(t 00 ) by repeated integration by parts. For the inhomogenous term, 
we have: 

llff-illi^ < ||<P-i||l2 

by definition of the multiplier, and since the Fourier support is contained in {|£| < 1/2}, 
we also have: 

Hs'-iIIlp < II^-iIIl 2 

for any p > 2 by Bernstein's inequality [2j Lemma 2.1]. 

Moreover, each of the (p^ is in S, and thus we have the convolution representation: 



<Pk(. x ) 



1 



7T 



9k(t) 



Thus we can write 



{ x -ty + {2- k y 

^2 P 2 -k * 9k, 
k>-l 



dt. 
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where the sum is convergent in H s if ip 6 H s . If ip 6 C a , then we have 

k 

<Y,2f> k \\<p k \\ L «, 

k 

< ||^|| c .^2^- Q ) fc . 

k 

We insert this expression in the variance: 



Var(Af n [<p]) = E(( ^ P 2 - fc * g k - E ^ P 2 _ fc * £ P 2 -< * 9l - E £ P 2 - ( * 
fc>-i fc>-i l>-i i>-i 

= Cov(7V n [P 2 - fc * g k ] Af n [P 2 -i * gi]) 

k,l>-l 

J2 f j 9k(t)gi(s)Cov(P 2 - k {- - t) P 2 -i (• - s)) dids, 



fc,«>-l 

which is formula (|10[). Thus, we may bound the variance as: 



11) Var(7VnM) <Y. \9k\\9i\Cow{M n [P 2 - k {t - ■)),Af n [P 2 -i(s - -)})dtds. 



k,i 

The following estimate will be useful below to control contributions from very hig 
frequencies. Its proof is a straightforward application of Theorem [7| 

Lemma 3.1. Let M > 1. Define the low-frequency cut-off operator Sm<P by 

Sm<P = <Pk- 
-l<k<M-l 

Then there is a constant C a such that, for <p £ C a , 

- SmMloo <c a 2- Ma y\\ ca . 

4. H 1+ CLT FOR GENERAL WlGNER MATRICES 

The goal of this section is to prove the following bound: 
Proposition 4.1. Let 

H := -^=W 

In 



be a Wigner matrix such that the entries of H satisfy condition (CO). Let 

z = E + in G T L , 

where 

T L := {z = E + in\\E\ < 5,0 < n < 10}. 
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We have 

(12) Var(TrG(z)) < ~ 

where G = (M — z)^ 1 is the resolvent for M. 

To prove the proposition, we will make heavy use of the following theorem from [T3] 
regarding uniform rigidity of eigenvalues. Before stating it, we introduce two definitions: 
first, we let m sc {z) be the Stieltjes transform of the semicircle distribution on [—2, 2]: 

2 /*2 

msc(z) = — \A — x 2 dx. 

7T J_2 Z - X 

We denote by m(z) the normalized trace of the resolvent: 

m(z) = — ti(H — z) -1 . 
n 

Theorem 8 (L. Erdos, E. Yau, J. Yin). Let H = (hij) be a hermitian or symmetric nx n 
random matrix, n>3, with Ehij = 0, 1 < i,j < n, and assume that variances are given by 
afj = E|/ijj| 2 = i. Suppose that the distributions of the matrix elements satisfy condition 
CO. There exist constants Aq > l,(f> < 1 such that for all L > Aologlogn and sufficiently 
large n, 

(13) P (u 2eSL {|m(z) - m sc (z)\ > ^fi ) < exp (-c(logn)^ 
where 

S L ■= {z = E + ir] : \E\ < 5, n" 1 (\ogn) l{)L < rj < 10}. 
For individual elements on the diagonal, we have, letting = max^ \ Gkk — msd, 

P(U 2SSL {A d > (l0S ^ 4L }) < exp f-c(logn)^ 



lnr\ 



Proof of Proposition 4-1 ■ Directly from Theorem [8], we can already conclude (12) for the 
region 

{z = E + irj: \E\ < 5, n~ 1+s < rj < n" 1/10 }, 

for any fixed < 5 < 1. We will now consider the regions < r/ < n~ 1+s and n -1 / 10 < 
77 < 10 separately, with 

o<5< 1 ' 



22 + e 

For the region r\ < n~ 1+l5 , we will use the fact that the eigenvalues are close to their 
classical locations. Let E = 3tz, and write 

n 

(14) TrG(z) = J2-, =^-- 

f-' Xi - E - ir] 

Consider a dyadic decomposition around E; define 

Uv = {j:\\ j -E\<n~ l+a }, 
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and for p = 1, . . . , log(n 1 a ) 

U p = {j : 2 p - 1 n~ 1+a < \\j -E\< 2 p n~ 1+Q } 

and 

Uoc = {j ■ \Xj -E\> 1}. 
a > will be chosen below. Given a sequence A n of events, we say that A n holds with 
overwhelming probability if 

P(A c n )<C k n- k 

for any k. By the local semicircle law, we know that for any 5' = e/100, we have \U P \ < 
2 p n a+s with overwhelming probability. In particular 

(E\U p \ 2 ) 1 / 2 <2^ +s ', 

since \U P \ < n. For eigenvalues with index in U p , we have the bound 

1 r < 2 1 - p n 1 - a . 

\*j-E\ ^ 

In Uq, we use the trivial bound 



1 i 

< rf ■ 



2 



\Xj — E — irj\ 

Split the sum (14) according to the partition {Up} to find 
E|TrG(z)| 2 < ri 2 - 2Q 2- fc 2~'E|[/ fc ||C/ z | 

+ — y^^nuowuki + \b\u 

Tj * — ' T) A 

Applying Cauchy-Schwarz in every term, we can bound this, up to a factor, by (n) 2+2<5 '(logn) 2 + 
n 2a+2S ' /rj 2 . Since we are assuming i] < n~ 1+<5 , E|Tr G(z)\ 2 is thus bounded up to a constant 
factor by r]~ 2 ~ e , provided we choose 6' < e/4, and a such that 

( + £ /4). (1 + ^)4 



We turn to the case where n -1 / 10 < r] < 10. In this case, we cannot obtain (12) 
from a direct application of Theorem [8j we cannot avoid the log factor appearing in 
large deviation-type bounds. To circumvent the problem, we estimate the variance of the 
resolvent trace using the set of self-consistent equations. 

We first introduce some definitions found in |14|. Let 



G ir .= {H-z)-\i,j), l<i,j<n. 

Define to be the (n — 1) x {n — 1) Wigner matrix obtained from H by removing the 
ith. row and column. We set 

G« = (ffW-z)- 1 (< >< 7). 
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The vector a 1 = (al)i<fc< n is obtained from the iih column of i/w by deleting its ith 
element. We denote by Ej the expectation with respect to the ith column of H, i.e. with 
respect to the random variables {hji)i<j< n . 

Lemma 4.2. The diagonal resolvent matrix elements Gu satisfy the following system of 
self- consistent equations 

Gu = 1 

-z-J2j a ijGjj + Tj 

where 

Tj := j4j + ha Z{ 

and 

(15) A = o-lG ii + J24^ 

Zi = zf- E^f, ZW = a < • G«a* = £ 4 ■ cfaj. 

In view of the previous lemma, we can write 

._ r _ 1 

-z - m sc - Q^. afjVj - Ti) 

Using the equation 

(m sc + z) = -m" 1 , 
and recalling afj = i, we can expand as 



4 E ^ - T *) - E ^ - T *) 2 + °d E u i - T *) 3 - 



r, - m scV 

n z — ' " n * — ' " n 

j 3 J 



We sum over the indices i and divide by n. Using the notation [v] := - Y^i v ii we obtain 
M = mUz)[v) - ^ Ti _ - T 4 ) 2 + 0((A d + T) 3 ). 

i i 

In the above equation, we have denoted 

T = max | T j | . 

i 

From now on, we will omit the dependency on z. The equation reduces to 

CM = C 2 -Vl 1 + C 3 - V(H - Tj) 2 + 0((A d + T) 3 ), 
n ^-^ n ^-^ 

i i 

where 



C\(z) x y/n + rj 
k = \2-Uz\, 

see jHj, (3.13), and C2, C3 are bounded independently of n. 
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are deterministic constants independent of n. We would like to prove the following 
bound: 



1 



Taking squares, then expectations, and using Cauchy-Schwarz, we have 

(16) E(|M| 2 ) < -L_ f 1e| E T *| 2 + i E KE(H - T ^) 2 ! 2 + E l°(^ + T ) 3 I 2 ) • 

We begin with the first and second terms in (16). The remainder term will be dealt with 
at the end. From Tj = Ai + ha — Zi, we have 

i iii 

The right hand side is bounded up to a constant factor by 



^(lX>l 2 + |V>MEz.l 2 ) 

\ i i i J 



Thus, it will suffice to establish the three estimates 

E|J> 



and 



E|J> 



E \J2 Zi 



|2 < « + r? 

i2 < « + V 
r*j 9 ? 



2 < W + 77 



The first one is obvious, since E|/ijj| 2 = -. We proceed to prove the second bound. Using 
the statement concerning individual elements of the resolvent in Theorem |8l we have 



(17) 



n 4^ 



GijG ji 



Gii 



■I J 



< 



m sc (z) 



n 



'J 



nJnrj 



For the first term, we write out 



u m (i)u m (j) 

Am £ 



where {u m {i))\<i< n is the eigenvector of H corresponding to the eigenvalue X m . We expand 
the double sum using orthogonality of the u m , and find that it is bounded by 



m sc (z) 



1 



// (A m - z) 2 



m sc (z)—m(z) 
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By the strong local semicircle law and Cauchy's integral formula, this is bounded by 1/ y/rj, 
up to a constant factor for rj > re -1 / 10 , with overwhelming probability. Squaring this, we 
obtain 

2 



m sc (z) 



n 



< I < K + 11 



V 



holding with overwhelming probability. Since \ Ai\ has a polynomial upper bound, we have 
obtained the desired estimate for the main term in (17). For the error term, we use 



\Gij\ < 



1 



nrj 



to obtain the bound, again using the fact that r\ > n 1 / 10 . 
It remains to bound the quantity 

e|£2.| 2 . 

i 

Since EZj = 0, the variance of [Z] is given by 

i n 1-1 

(18) - 2 B\ = z« z e + ^ E E \ z ° 



n- 



11- 



We shall focus on the case where a = 1 and (3 = 2. Writing g[V as 



G 



(i) 

kl 



p(12) p (l) 
~r r kl ■ 



where P^ and P^ (for k,l^l) are defined as 



,(12) 



(12) p (l) ._ ^fc2 ^21 



• — ^kl ' ^feZ 



G 



(1) 

22 



if k,l^2. 



P 



(12) 



0,P 



(1) .. 



G 



(i) 



if k = 2 or I = 2. 



With these notations and denoting IEj := I — E a . , we can express 

Z 1 = IEia 1 • P (12) a* + IEia 1 • PV 

and 

Z 2 = IE 2 a 2 • P (21) a 2 + IE 2 a 2 • P 2 a 2 . 



Using the previous expressions, we see that to estimate the right side of (18), it suffices 
to obtain the following bound: 



Ela^pWa 1 ! 2 ^ 



1 



(nrj) 2 
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Follow the idea in section 8.1 of [15], we write 
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Ela^pWa 1 ! 



E 



E4 



a l 



k.l 



G 



(i) 

22 



Expanding, we have 



\ " ~T 1 \-^Gk 2 G 2 \G v2 G 2q 

EE i E a l a i a v a \ 



k,l,p,q 



\G 22 \ 2 



There are only a few cases where the first expectation is nonzero. One such case is k = 
I, p = q. The other situations can be dealt with in similar fashion. Taking the first 
expectation in, we have to bound 

/ \ 2 



■n- 



\G- 



22 



We expand 



IG22I = m sc (z) + O 



(logn 



AL 



nrj 



with overwhelming probability. Thus, the only important quantity is 

2 

(19) 1 



^2 E (j2 G k2 G 2! )j 



Expanding the above sum, we obtain 



n 2 E ( E 



:E 



dz 



G 22 (z] 



Using Theorem [14J and Cauchy's estimate once more, we obtain 



n 2 E ( E 



Gk 2 G 2 k 
~]G> 



22 



< J_ 



with overwhelming probability. Again, all quantities involved have polynomial bounds, so 
we have obtained the desired estimate. 

Next, we direct our attention to the second term on the right hand side of equation (16). 
Using Theorem (|8j), we see that the term ^(X]Jw] 2 ) 2 will be bounded by n -4+e ?7 -4 and 
since r/ > n -1 / 10 , will in turn be bounded by n _2 r/ _1 if rj > n~ 2 / 3 . So it remains to bound 
the term 

E ^i E T *i 2 ^ ( E i E ^ 2 + E i E ^i 2 + E i E ^ 2 i 2 ) • 

i \ i i if 



IS 
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Again by Theorem ([8j) and equation (15), the first term on the right hand side is harmless. 
The second term is trivial to bound so we focus on the third term. Applying Cauchy- 
Schwarz twice, we have 

2 



(20) 



1 



The desired bound then follows immediately from lemma 4.3 below, the definition of Z, 
and the estimates 

1/2 



£k#n < 



n. 



By (3.39), (3.41) in [Hj, we have 

1/2 



W|2 

ij 



< 



TTTfe) 



1/2 



x (fe) _ . 



holding with overwhelming probability. The proof of the next lemma can be found in 
Appendix B of [16] . 

Lemma 4.3. Let ai (1 < i < n) be n independent random complex variables with mean 
zero, variance a 2 and having uniform subexponential decay. Let Bij E C (1 < i < n). Then 
we have that 

n n n 

F{\ J2 <kBuai- <r 2 B u \ > Da 2 (Y^\Bu\ 2 ) 1/2 } < Cexp(-cD^) 



8=1 



1=1 



and 



J2aiB ijaj \ > Da 2 (J2\ B ij\ 2 ) 1/2 } < Cexp(-cD^+ s i 



Lastly, the last term in equation (16) can be bounded by n - 3 / 2 + e ^- 3 / 2 because both A^ 
and T are bounded by n~ 1 ^ 2+e r]~ 1 ^ 2 with overwhelming probability by Theore m M For 
rj > n -1 / 10 , we have ra _3 / 2 + e ^- 3 / 2 < n~ 1//2 ry -1 / 2 and this concludes the proof of (|12T □ 

We are now ready to prove Theorem [lj Write ip = XV + (1 — x)^P^ where x is a cut-off 
function equal to 1 in [—3,3] and equal to in [— 4,4] c . Then we have 

Var(AA n [99]) = Var(M n [ X <p] + Af„[(l - x)<p]) 

< 2Var(AA n [ X ^]) + 2Var(AA n [(l - x )<p])- 

We bound the second term by the second moment E|A/" n [(l — x) ( / J | 2 > which can be estimated 
as follows: 

P(l E (C 1 " Ml > 0) < P(Ai > 4) + P(A n < -4). 



REGULARITY CONDITIONS IN THE CLT 19 

The last quantity is bounded by e~ nC (|16j. Lemma 7.2). Thus 

EK[(1-^]| 2 = 2 / xP(| Y>(Ai)| >x)dx 

i 

<n 2 [Mli«-P(|£p(AOI>o) 

i 

2 — n c II II 2 <- II ||2 

= n e \m\ L oo < \\<p\\ H i+e- 

The final inequality follows from the Sobolev inequality. For the first term above, use the 
decomposition (fTTj) for if, but not for X- 



Var(Af n [ X <p]) <Y J j l<7fc||Si|Cov(Af n [x(0^2-*(* " -)],^*n[x(-)^(« " OD^ds. 

Apply Cauchy-Schwarz to the covariance on the right side of the previous equation. Since 
the integrals are now decoupled, we concentrate on bounding: 

Y J |^|Var(AA n [ X (-)P 2 ^(t--)]) 1/2 dt. 

k 

Apply Cauchy-Schwarz again to the integral 

(21) 1 2 {1+e)fc b fc |2-( 1+ ^Var(AA„[ X (-)P 2 - fe (t - -)}) 1/2 dt 

k 

and obtain 

Y (/ 2( 2 + 2 ^| 9fc | 2 dt) 1/2 (| 2-( 2 + 2 ^Var(AA n [ X (.)P 2 - fc (t- -)])dt) V2 . 
We are left to examine the quantity: 



/ Var(A/- n [x(-)P 2 -*(*--)])d*- 



Split the integral: 

(22) f Var(Af n \x(-)P 2 -K(t -•)])&+ [ V a r(Af n [ X (-)P 2 - k (*-•)]) dt. 

The integrand in the second term is the variance of a smooth function with support in 
[—3, 3], and such that 

C 

||X-P 2 ~ ft (* - -)llc 2 ([-5,5]) ~ _ 

By M. Shcherbina's variance bound [29] for linear statistics of H 3 / 2 functions, this term is 
thus integrable. The first term in {1221) can be written as a sum: 



Var(A/k[x(0P2-*(*-0];N,|A n | < 2) 

+ Var(AA n [ x (-)P 2 - fc (t - •)]; {N, |A n | < 2} c 
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Since x = 1 on [—2, 2]: 

V a r(Af n [ X (-)P 2 - k (t - •)]; |Ai|, |A„| < 2) < Var(AA„,[P 2 - ft (t - •)]), 



which is bounded by 2 2k+e by (12). On other hand, 

Var(AA„[ X (-) J P 2 - fc (i - ■)]; {N, \K\ < 2} c ) < 2 2k n 2 e~ nC 



Thus, (21) is bounded up to a constant factor by: 

1/2 



Choosing e suitably, a final application of Cauchy-Schwarz and Theorem [7] end the proof 
of the bound ([6]). The central limit theorem then follows from the discussion in Section [2] 

Remark 1. The same proof can yield the CLT for functions in for all p > 1 in this 

case. 

5. Variance bound and CLT for C^ + : GUE case 

Our starting point is the following representation for the variance, based on the deter- 
minant kernel representation of the 2-point function: 

(23) Var (Af n [tp]) = - ft \<p{x) - ^{y)\ 2 (K n (x, y)) 2 dy&x, 

with the kernel K n being defined in terms of Hermite functions by the Christoffel-Darboux 
formula: 

Kn(x,y) = 2^^3K x )VAy) = • 

i=o x y 

The Hermite functions are defined by: 

fl*(x) = (-l) fc e- 8 ^ 
j.( x ) - II M e - 2 /2 

m } ~ (27T)V4(fc!)l/2 

, . , /n\ V 4 ~ / n 1 / 2 x\ 

m*) = y — ■ 



Let us recall how formula (23) is derived. We can write the variance as: 

Var {NrM) = EAT 2 [p] - (EW„M) 2 

= E^^(A i ) ¥ ,(A j ) + E^^ (A . ) 

i^j i 

-^E^(Ai)E^(A,-). 



hi 



REGULARITY CONDITIONS IN THE CLT 21 

In the second equality we have expanded the definition of M n [<p\ . By definition 

(24) EJ2<f(X i MX j )=n(n-l) fj tp(X)ip(p) pffl (A, p) dXdp. 

(25) E^(A,) 2 = n^ 2 (A)p«(A)dA 

p 2 {X) W(A)p«(/i) - n(n - l)p^ 2 \X,p))dpdX 
= ^ // (vAA) + p 2 {p)) (n 2 p«(A>«(/i) - n(n - l)p^ (X, p))dpdX, 

(26) ^E^AOE^A,-) = n 2 || tp^X) (X)pW (jj) dpdX. 

(k) 

In the above, pfi denotes the fc-point correlation function of the GUE eigenvalues. By the 
Dyson-Gaudin-Mehta computation, we have: 

p£\X)=K n (X,X) 
p {2 \X,p) = K n (X,X)K n (p,p) ~ K n (X,p)K n (p,X) 



Using symmetry of K n , i.e. K n (X, p) = K n (p,X), we arrive at the variance formula (23). 

Let us return to the variance bound. M. Shcherbina and L. Pastur [28], Theorem 5.2.7, 
show that 

(27) Var^MHjL/Y M^fag ^- ^d^ + offl 

for Lipschitz functions ip. By carrying out the computation carefully, we will see that the 
result extends to C 1 ^ 2+ functions, and the convergence is uniform over bounded subsets of 

Let us first show that the main term in (27) can be controlled by \\p\\ i +e . The first 



step is to make use of the C l l 2+t hypothesis. Let e > such that p 6 C 1 ^ 2+e , then 
(p(x)-p(y)) 2 A-xy 



-2,2]2 ( x ~ V) 2 \fl - x 2 ^/A - y 



dxdy 



Thus we just need to check that the integral is finite. It is clear that the integrand is locally 
integrable away from the points (—2,-2) and (2,2). Also, the integrand is positive. By 
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symmetry, it is sufficient to check finiteness of the integral: 
1 4 — xy 



o Jo Ix-yl 1 ' 2 ' V4-xV4-y 2 



dydx 



1 f 2 f x 1 A-xy 



~2j J \ X -y\l-2e^2^^2^ dydX 

Writing 4 — xy = 2(2 — x) + (2 — the last integral can be bounded by: 

Since \/2 — y > \/2 — x on the domain of integration, the first integral is convergent; by 
considering the cases \x — y\ < |2 — x\% and \x — y| > |2 — x\*, we see that the second 
integral is convergent as well. 



asymptotics for the Hermite polynomials 



We now turn to the computation of (27). The main tool will be the Plancherel-Rotach 

[35] . Chapter 8. We begin by showing that we 



can replace the integral in (23) by an integral over the square 

max(|x|, \y\) < 2 + n~ v 
for some < v < 2/3 to be selected later. This follows by the simple bound [28] (5.2.8) 

K n {2 + sn- 2 '\ 2 + sn~ 2 / 3 ) < C^W 08 "". 

This bound in turn follows from Plancherel-Rotach asymptotics in the transition region. 
Applying it, we find that 

Var (Af n [(p\) = \ II \<p{x) ~ V(y)\ 2 (K n (x, y)) 2 dydx + M^OipT 00 ). 

J J\x\,\y\<2+n-" 

We will now further replace the integral over the square {|x|, |y| < 2 + n~ u } by an integral 
over {|x|,|y| < 2 — n~ u }. The difference between these two squares is contained in the 
union of four rectangles. Without loss of generality, we can consider the integral over the 
"bottom" rectangle 

R n = [-2 - n~ u , 2 + n~ u ) x [-2 - n~", -2 + n~ v ]. 

We write: 

(p(x) - ip(y)) 2 



, N2 -(^n{x)^n-l(y) ~ ^n(y)tpn-l(x)) dxdy 

(x - y) 1 



/ I nz^(ipn(x)ip n -i(y) - ip n (y)ipn-i(x)) 2 dxdy. 

-2—vr v J—2—n~ v \ x ~ V\ 

The last integral is bounded up to a constant factor by: 



2+n~ v r -2+n~ v 



1 



j^2e ((lpn(x)lp n -l(y)) + {^n{v)^n-x{x)) ) dxdy. 
2-n _ " J-2-n~ v \ x ~ V\ 
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By symmetry (because i\) n and ip n -l have similar behaviour as far as we are concerned), it 
suffices to deal with: 



(28) 



2+n - " 



2-n~ v 



-2+n,-" 



1 



2-n-" 



\x-y 



l-2e 



dydx. 



Using the bound < Cn 12 for all n [35], Theorem 8.22.29 we have (VVi-i(y)) 2 < 
£< n -i/6 n i/2 _ (j n l/3 l nser ting this bound in the inner integral and computing, we find: 

-2+n-" 



(29) 



(V„-i(y)) 2 , 

2-n"" \X — V 



rz^dy < C e n 1/3 — —— T 3^1[_ 2 _ 2n -- _ 2 +2n--]=(^) 



n 

\x + 2 1 

+ C7 e n 1/3 l [ _ 2 _ 2n -, i _ 2+2n -,i(x)n- 2 '"" 



-2-2n"", -2+271-"]' 

The last estimate is obtained by direct computation. For the first term on the right, 
we have used the Taylor expansion of (2 + n~ u /a) 2e with a = x + 2, valid in the range 
rT v < ||a|. We insert this bound in the x integral in (28). In the region —2 + 0(n~ u ), 
corresponding on the right side of (29), we use the bound (ip n (x)) 2 < Cn 1 / 3 , leading to a 
bound of order 

C t n 2 l 3 n~»~ 2 ^. 

Choosing, as we may, 

2 I 2 

< v < -, 

31 + 2e 3 

the contribution from the second term to (28) is 0(n~ c ) for c > 0. The first term on the 
right of (29) leads to a contribution bounded by: 



(30) n 1 ^n- v I (Vn(x)) 2 |x + 2|- 1+2e dx 

-2+2n~" 



,1/3-" 




2+rt - 



For the first integral, notice that since v < 2/3, we are far clear from the transition region 
of the rescaled Hermite functions i/j n , and so we can apply the asymptotics in [IB] (see also 
(5.1.9)), which show that ip n , tp n -i are bounded by 

C 



v / 2^2 

in the region of integration. Thus the first term on the right side of (p0| is bounded by: 



1 1 



' 2+2n- V / 2^k + 2| 1 - 2e V ' 



For the second integral on the right side of (30), notice that \x + 2| is bounded below on 
the range of integration, and that on the other hand, the total integral of (ip n (x)) 2 is equal 
to 1, Thus the contribution from this integral is 0(n 1//3_y ). 
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We have so far established the asymptotic behaviour: 

War (M n [<p]) = // \<p(x) -<p{y)\ 2 (K n {x,y)) 2 dydx + C\\ip\\ 2 cl/2+t n- c . 

J J\x\,\y\<2-n-" 

We now wish to replace n~ v (whose smallness has been useful so far) in the integral range 
by a fixed 5 > 0. We will bound the difference, which is the sum of four integrals of the 
form 

/ / 

J2-8<\x\<2-n~" J\y\<2-n" 

in terms of 5, independently of n. Indeed, as explained above, since v < 2/3, the Hermite 
functions are bounded in the region of interest. Inserting these L°° bounds in the integral, 
and extending the upper limits of integration from 2 — n~ v to 2, we find that the difference 
is estimated by 

(31) CM^ 1/2+E / / 1 dxdy=[cp] 2 cl/2+€ -0(6). 

J2-S<\x\<2 J\y\<2 \ x ~ U\ 
Thus we are now concerned with computing the quantity: 

(32) // \<p{x) - <p(y)\ 2 (K n (x,y)) 2 dydx 

■>J\*\,\y\<2-6 



for 5 > 0. We need to show that it matches the first term in (27) with an error uniform in 
||(^|| c ,i/ 2+e < K for any K > 0. Expanding (K n ) 2 , we find: 

(lpn{x)) 2 {lpn-l{y)) 2 + {^n{y)) 2 ^n~l{x)) 2 ~ 2tp n (x)lp n (y)^n-l (x)lpn-l (y) ■ 

By symmetry, it suffices to consider the contributions of the first and last terms. Using 
the Plancherel-Rotach asymptotics and trigonometric identities, we have: 

(M*)) 2 = y=p( cos Va(#) " 9/2 - vr/4) + O^" 1 )) 

' (l + cos(2na(^)-0-7r/2)) + O(n~ 1 ), 



uniformly in |x| < 2 — 5, where 

x = cos 6 
a(6) = 9-sm(29)/2 

Similarly 

ipnix^n^ix) = 1 (cos(2fl) + cos(2na(^) -9- tt/2)) + 0(n~ l ). 



V2 



X 



We will now simply insert these expansions into the double integral over |x|, |y| < 2 — 5, 
and use the C a smoothness of the integrand (as opposed to the Riemann-Lebesgue lemma 
for L 1 as in Pastur and Shcherbina's computation) to obtain a uniform rate of decay (for 
each 5 > and uniform in Hvllcn/a+e < K). 
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We start with the integrals resulting from the square terms (tp n (x)) 2 (tp n -i(x)) 2 and 

(My)) 2 ^n-i(y)) 2 : 

L^ (*fo)-y (1 + C os(2na(#(x)) - 6{x) - vr/2)) 

x (1 + cos(2na(d(y)) - 9{y) - vr/2)) dxdy + ^^Offn- 1 ). 

The Os(-) term is justified by | cos(-)| < 1 and by substituting \x — y\~ 1+2e for the ratio 
{<p{x) — tp(y)) 2 /(x — y) 2 . Similarly, the integrals resulting from expanding 

- 2 tp n (x) ip n (y ) 1p n - 1 (x) ip n (y ) 

gives a contribution of 

1 (y(x)-<^)) 2 + co<2na{e{x)) _ e{x) _ v/2)) 

'{|ar|,|v|<2-5} V4 - X 2 a/ 4 ~ V \ x ~ V) 

x (y + cos(2raa(#(y)) — #(?/) — vr/2)) dxdy + [y]ci/2+e 
Adding the main term in the last two displayed equations, we find 

4-xy (ip(x) - ip{y)) 2 



{\x\,\y\<2-5} V4- x 2 y / 4-y 2 {x-yf 



-dxdy, 



as expected. It remains to show that the oscillatory integrals decay uniformly for ip in a 
bounded set of C l l 2+t and each 5 > 0. For this, it will suffice to deal with integrals of the 
form 

( oo, f 2 ~ 5 f0 (a < xy>(*)-P(v)) a dx 

(33) / cos(2na(#(x)))- 



2+5 ' "' (x-y) 2 \fl 



X 



2 



a and x are related by the change of variables 

^M = -2sin0(x) = -(4-x 2 ) 1/2 , 
dx 

non singular away from x = and x = ±2. We treat ( |33[ ) as an oscillatory integral with 
non-stationary phase and Holder continuous amplitude. The guiding model is the integral 

(34) J(t)= / e**^*A(x)dx 



with |(/>'| > e and an amplitude A £ C a . This can be dealt with by writing: 

i<t>{x)t a l x+ ZL^ )dx= / e^V*" A(x) dx + 0|| a |, £ (l/t) 

V t^(s)y Ja+0(l/t) 
rb 



J e i(tKx)^) e iO*{i/#) A { x } dx + o Moo>e (l/t) 

-l(t) + 0(l/t). 



26 PHILIPPE SOSOE AND PERCY WONG 

Thus we can write 



2|/(t)| 



b 



M.x)t (A(x) - A(x + ir/(t(f>'(x)))) dx 



+ 0(l/t). 



By the a-H61der condition, the main term in the last equation is bounded by C||a||cc* -t~ a if 
a, b < oo. A similar analysis around around points x where 4>'(x) = but <j)"(x) > e shows 
that the contribution to the integral from a neighbourhood of that point is C||a||cai _a / 2 . 

We will apply the same strategy to ( |33| ), but we have to deal with the fact that the 
amplitude factor (seen as a function of x) becomes singular in a neighbourhood of x = y. 
Away from y, the algebra property of C a ensures that the C a norm of 

(^)-^)) 2 1 



(x - y) 2 yfl~ 



x* 

has a bound depending on 5, [<p] c i/2+e and the distance to y. Let us roughly calculate the 
dependence on the distance to y, i.e., we are interested in giving an estimate for the local 
norm 

II " HcV2+e(i) = 115/ • llc ,1 /2+6, 
where / C (—2 + 5,2 — 5) and dist(I, y) = s, g is a C°° function supported on 

{dist(2» < s/10} 

and g = 1 on /. The function x i— > (f(x) — ip(y)) 2 inherits its Holder norm from (p, since: 

\(<p(x) - v{y)f - (v(s') - v(y)f\ < \<p(x) - <p(x?)\ ■ MMoo- 

We may ignore the factor l/vi — x 2 , since its C 1 norm is bounded for |x| < 2 — 5. To 
estimate the C a norm of — y) 2 , we simply interpolate between the L°° and Lipschitz 
norms. For x,y £ I + -6(0, s): 

^M<i/(x)-/wi-^Mr 

\x — y\ a \x — y\ a 

<c\\f\\ l - a \\f\\ a ca . 

If dist(J, x) > s, we thus have 

\\l/(x - y) 2 ||Loo (/) < s -2 , 

W(z - yfWcni) Z - y| 3 llL-(/) < s~ 3 , 
l|i/(^-y) 2 lbv 2+E <s- 5/2+e . 

On the other hand, away from x = ±2, only the squared factor is singular, and so we can 
bound the integral by 

CiMcV2+ £ / I 1 |1 _ 2e dx < C s s 2e . 

J\x-y\<s,\x\<2 \ x ~ V\ 



It is now clear how to proceed: we split the integral (33) according to the distance to y: 

r 

J \x\<2-8,\x-y\<n~ a J \x\<2-5,\x-y\>n- a 
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with a > a parameter to be adjusted later. The contribution from the first term is 
bounded by C5[(^]^ 1/2+f n _a2e , as just explained. For the second part, we can (almost) 
repeat the computations leading to the estimate for the model integral (34), using the 
estimate Cs[<p]Qi/2+ c n a ( 5 / 2 ~ eS) for the 1/2 + e Holder norm of the amplitude. All in all, we 
are lead to an estimate of the order 

However, we have ignored the fact that our L°° norm estimates do not allow us to bound 
the error at the endpoints resulting from the translation of the variable x, i.e. the integral 
over 

[y-n- a - Cn-\y-n- a }. 
If instead of using the L°° bound, we replace the integrand by \x — y\~ 1+2e once again, 
we obtain an error of order M^i/2+ e n_2e - Choosing a appropriately (a = ^§^r), we find 
that the model oscillatory integral ( 33 ) decays at a power rate uniformly for <p in bounded 
subsets of C l l 2+€ and uniformly in \y\ < 2 — 5. We can use the scheme above to deal with 
all the oscillatory terms we encounter, and thus we have established that 

Var(A4M) = // ~ ^ dxdy 

+ [^ /2+ MS) + Mcl/2+e>6 (n- K ) 

for some k > 0. 

Given K such that [|y||c?i/2+e < K, we can choose 8{K) to make the error in (31 ) smaller 
than e', and then choose n(5) so as to make (32 ) e'-close to the main term in the asymptotic 
(27). Thus for n > n(K), we have 

which is enough to close the approximation argument leading to the CLT. 



6. CLT FOR C l / 2+ : A CLASS OF WlGNER MATRICES 

In this section, we prove the central limit theorem for C 1//2+<E test functions for random 
matrices of Gaussian convolution type. To achieve this, we will use a comparison procedure 
for the resolvent as in [33]. Let us first precisely describe the matrices we will be working 
with. The matrices we will consider have the form 

(35) V^M = \{W + V), 

where V is an independent GUE matrix and the distribution of the entries of W = 
( w ij)i<i,j<n satisfies condition CO (see ([5])). The entries $lwij, Sswij, 1 < i < j < n 
are all independent and the first five moments of the complex entries Wij match those of 
the GUE. That is 

E[($lwij) a (5Swij)P] = EX^I 3 
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for any choice of < a, (3 < 5 such that a + /3 = 5, where X and Y are independent iV(0, ^) 
random variables. We remark that the entries of ^/nM have sub-exponential tails. 

We denote by (fk the Littlewood-Paley projection of ip onto frequencies of order 2 k (see 
Section |3j). Write: 

v = E + E + E 

\£:<logn log n<fc<2 log n k>2\ogn J 

If G C 1 / 2+e , we use the trivial bound 



and Lemma 3.1, and obtain 

EKV' n ]| 2 <^n- e . 

Thus, the sequence A/^fy? 3 '™] is a remainder whose second moment tends to zero. The 
contribution to A/"°[i^] from ip 2,n cannot be estimated quite so crudely. Nevertheless, we 
will show below that 

(36) E|AH^ 2 >1 2 < n" ce 

for some c > 0. Thus (p 2,n does not contribute to the asymptotic distribution. 

Assuming (36) for the moment, it suffices to show that Nffif \ ls asymptotically normal, 
with the correct limiting variance. For this, we will compare A^ft/? 1 '™] to the corresponding 
quantity defined in terms of GUE eigenvalues. Let 

n n 

3=1 3=1 

where Xj, 1 < j < n denote the eigenvalues of a GUE matrix normalized by l/\/2n. We 
express the variance in terms of the GUE quantity: 

Var(A/;V n ]) = Var(A/;° GUE [^"]) + (Var£A£ gue^^D " Var(AA°[ v? 1 -™])) . 
By Q, the first term is bounded. Using equation (10), we estimate the second term by 

(37) n 2 Y, JJ\C(ztzi-,M)-C(zlz l 2 -,GVB)\\g k (s)\\ gi (t)\dtds. 

— l<fc,i<log n 

Here, z\ = s + i • 2 _fc , z\ = t + i- 2~ l , and 

C(zx,z 2 ;M) = Cov(3a(*i), Zs(z 2 )), 
C(z 1 ,z 2 ; GUE) = Cov(9sgue(^i),9sgue(2 2 )), 

s( z ) = -tr(M - zY 1 , 
n 

sgve(z) = - tr \-^=V - 
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We split the integral in (37) into two regions: 

(38) n 2 [[ \C(zf,zi ] M)-C(zt4-,GVE)\\g k (s)\\g l (t)\dtds 

J J [—3,3] x [—3,3] 

+ n 2 [[ \C(ztz 2] M)-C(z k 1 ,z l 2 -,GVE)\\g k (s)\\g l (t)\dtds. 

J J ([— 3,3] x [— 3,3]) c 
We first show that the integrals 

(39) n 2 [[ \C(z k ,z 2] M)-C(z k ,z 2 -,G\JE)\\g k ( S )\\ gi (t)\dtds 

J J [—3,3] x [—3,3] 

are bounded by H^dcrx/a+e, with factors summable in k and I. 

To compare C(z k , z 2 ; M) and C(z\, z l 2 ; GUE), we shall make use of the following slight 
modification of a result of T. Tao and V. Vu 



Proposition 6.1 (Tao, Vu). Let Mq be a random Hermitian matrix where each entry 
satisfies condition (CO), and let £ be a random variable independent of Mq. Let 

M = Mq + —f=£A, 



where A is an elementary matrix. Denote the Stieltjes transform of the empirical eigenvalue 
densities of Mq and M by sq and s, respectively. Then, for fixed rn > 1 and uniformly over 
\E\ < 3,rj > n -1-5 for some sufficiently small 5 > and with overwhelming probability, we 
have 



(40) s{E + in) = s (E + irj) + Y J i 3 ^ j/2 c J {7 1 ) + o( n -("»+i)/2(i-C«)) min{l, — -}, 
where 



Cj (ri) = 0(ra c ^min{l,— }). 

nrj 

The proof of the above proposition follows as in |33j . To derive the main result in 
Tao and Vu only needed to compare resolvents for a fixed energy E. Using Theorem [8j 
together with the bound 

^ss(E + irjz) ~ rji ' 

Lemma 16 in [33] will hold uniformly over the range \E\ < 3. By |33j . Proposition 13, this 
yields the expansion above. In particular, the condition that rj > n~ 1-5 is crucial. 
Setting 

m 

S(m,rf) = ^n~^' 2 Cj{rj) 
i=i 
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and taking imaginary parts on both sides of equation ( 40 ) , we have 

(41) E(Zs(z k )Zs(z k )) = E(9too(*f)3so(4)) 

+ E3?s (^) • 9 (S(m, 2~ l ) + 0{n~^ {l+cs ^) ■ min(l, n~ l 2 l ) 

+ E9s (4) • 9 (S(m, 2~ k ) + 0(n~ m ^^ 1+c5 '>) • min(l, 



+ E9 ( 5(m, 2"') + 0(ra-^ (1+c<5) ) • min(l, n -1 ^ 



x 9 (S(m, 2- k ) + 0(n- 2 £ 1 < 1+CT )) • min(l, rTH k )\ . 



In the above, we have omitted error terms from exceptional sets whose complements have 
overwhelming probability. This is justified since all quantities are polynomially bounded in 
n for k, I < Clogn. We will apply the above estimate in the situation where variables £i, £2 
take the place of £ and the moments of £1,^2 match up to fifth order. When we take the 
difference, the terms involving moments of fifth order or less can be ignored. However, we 



cannot afford the error terms appearing on the second and third lines of (41 ). On the other 
hand, as will be seen below, the factor min(l,n — 1 2 ) • min(l,n _1 2') appearing in the final 



term on the right side of (41) is sufficient to allow control of the difference of covariances. 



We use (40) in the second term E(s(zi))E(s(z2)) to obtain the cancellation we need: 



(42) E(3s(sf))E(9to(4)) = E(3 ao (*f))E(3*o(4)) 

+ E9s (z 1 fc )E € 9 (s*(m, 2~ l ) + 0{n-^ l - cs ^) min(l, rT^ 1 ) 

+ EQ?s (4) E £ 9 ( s ( m , 2 ~ k ) + Oin-^+C 1 - 05 ^) min(l, n~ 1 2 k ) 

+ E f 9 (S(m, 2- k ) + 0(n- m ^ 1+ ^ 1+cs y) min(l ) n- 1 2*)) 

x E € 9 (S(m, 2~ l ) + 0(n-^ {1+c5) ) min(l, n~ l 2 1 )} . 

Because so does not depend on £, if we first take expectation with respect to £, we see 
that the terms on the second and third lines of ( 42 ) cancel the terms on the corresponding 
lines in (41) to any desired order provided m is chosen large enough. If we take m large 
enough, say m > 8, we can gain enough powers of n to make the error term on the second 
and third lines less than 



n 



-3+e' 



min(l, 2 k /n) min(l, 2 l /n), 



for some small e' depending on 5. We have now established that for two matrices M, M' 
differing in one pair of entries whose distributions have 5 matching moments and for some 
small e' > 0, 



(43) \C(z k ,z l 2 ;M)-C(z k ,z l 2 ;M')\ < 0(n~ 3+e ') min(l, n~ 1 2 k ) min(l, rT 1 ^ 
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We have to swap 0{n 2 ) entries to bridge between the GUE and our Wigner matrix. 
Applying equation (43) repeatedly and taking into account the factor of n 2 in equation 



(39), we see that if two matrices match up to 5 moments, then the difference in the 



covariance is bounded by 

(44) n 1+e ' ■ minj 1 , — — r ) mini 1 , — — r } • 



The estimate (44) allow us to close the argument. To see this, we sum over k,l of (39) 



and using rt44h , we have 



n 2 



V // \C(zlzi;M)-C(zlz l 2 ;GHE)\\g k ( S )\\ m (t)\dtds 

-l<M<lognM-3,3]x[-3,3] 



^2 \\9k\\L°° min(l,n l 2 



, — Kfc<logn 




Inside the bracket, we have 



H^IU-mina^rr 1 ) = £ \\g k \\ L oo2 k ^2 k ^) . rT 1 

-l<Ai<logn — l<fc<logn 

<n- 1 sup(2^^)|| ffi || i o ) 

3 

— l<j'<log n 

< n~2~ e ' \\(p\\ 01/2+e . 

Squaring both sides, we gain a factor of n~ 1 ~ 2t which allows us to cancel out the factor 



of n 1+e ' provided we pick e' sufficiently small in (43), thus ensuring that the error stays 



bounded up to a constant (in fact o(l)) factor by || • I n conclusion, we have bounded 
VarCA/^ 1 '™]) by C7||^" 2 



C-l/2+6 



We now deal with the integrals in (22) over the region ([—3,3] x [— 3,3]) c . We illustrate 



the procedure with the integral over the infinite rectangle 

{(s,t) : (s,t) e [-3,3] x [3,oo)|, 



the other cases being similar. We use the deterministic version of (40) in [33], as we have 



done in the case s,t € [—3,3] above. However, instead of bounding the coefficients Cj(rj) 



as in Proposition 6.1 we return to the expansion in [33], Proposition 13. We expand both 
( sss(z k ), ^ss(z l 2 ) to 8th order as above, and perform all the same cancellations, but use a 
different bound for the coefficients appearing in the expansion of s(z l 2 ), the resolvent at 
energy t G [3, oo). The estimate provided by [33] for Cj(rj) is: 

Cj = 0(\\Ro\\ J {ooA) mm(||flo|| (DOil ), 1/nrj)), 
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where Rq is the resolvent of the unperturbed matrix Mq in the statement of Proposition 



6.1 and 



ll-Ro||(oo,l) = m a 5 \(Ro)ij\ 

l<i,j<n 

On the event max{|Ai|, |A n |} < t/2, we have 

1 „ 1 



max 



< 



Wi - 1 ~ t 



for t £ [3,oo). By [16], Theorem 7.2, the complement of this event has probability 

£ e -n c logi < 3 -n c /2 t -n c /2^ 

The conjunction of the event and derealization of eigenvectors, a consequence of the local- 
semicircle law, implies, as in |33j . Lemma 16, that 

Cj{rj) = 0{n e min(n e t~ l ,1/nrj). 

uniformly in r] > n -1-5 . Using these observations, we obtain the bound: 

|C(zf, 4; Af) - C(^i> 4; < 0(n- 3+e ')t- 2 min(l, n -1 2 fc ) min(l, n -1 ^), 

for s,t£ [—3, 3] x [0, 3), which allows us to perform the integral over the infinite range and 
conclude the argument as previously. 

It only remains to justify equation ( 36 ) . Only here do we utilize the Gaussian convolution 
structure. We begin by stating the following results due to E. Brezin and S. Hikami [3], 
and K. Johansson |22j . (Proposition 1.1, Lemma 2.1, Lemma 3.2). 

Proposition 6.2 (Brezin-Hikami, Johansson). The symmetrized eigenvalue measure on 
R n induced by the Johansson matrix M defined in (Sty has a density p n (x), given by 



(45) Pn{x)= q aVn (x;z(H))dpW(H), 
where a 2 = \ and 

q a y n [x,z) An ^Y et ^ e h,k=v 

A n = Y\i<j I x i ~ x j I * s ^ e Vandermonde determinant. 

dp( N \H) denotes the measure on the space of Hermitian matrices induced by the entry 
distribution of the matrix H = -^j^W (see (35)), and z = {zj{H)}^ =l is the vector of 
eigenvalues of H. 

Proposition 6.3 (Brzin-Hikami, Johansson). The correlation functions for q a 2 i n {x\ z) are 
given by 

n! f 

(46) R n m (xi, . . . ,x m ;z) = ^7/ q a 2/ n (x; z)dx m+1 . . .dx n 

(n-m)\J R N-m 1 

(47) =det(Kf/ n (x i ,x j ;z))% =1 , 
for some kernel K n . 
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In equation (2.19) in [22J, Johansson gives an alternative form for the determinantal 
kernel appearing in the previous proposition: 

lC n (u, v; z) = e^^^-^Kt^^ v; z), 

where u) is any constant, IC n (u,v; z) also serves as a determinantal kernel which produces 
the same correlation functions q a 2j n (x; z). 

We use the following variant of the asymptotic analysis of the diagonal kernel Ki n {u, 2) 
in [21] , derived using ideas from [13] • The proof will be found in the next section. 

Lemma 6.4. Let 5 1 ,5 2 > 0. Let n R , n = {s G C : \Us\ < R,^ 1 ^ < |Gs| < R}. Let Z n be 

the set of z G M. n for which the local semicircle law holds, that is, the set of z such that the 



event on the left side of (13) is satisfied with m{z) = - z _ zi . Then, for z G Z n and 

, we ha? 



any x such that \x\ > y/2 — n 1 / 3 + 5 2 ; we have the estimate 

1 



r^C>l,2 " 



np sc (x) 

Here, p sc denotes the semicircle density: 

p sc (x) = -^2-x 2 l { _^^{x). 

Assuming the lemma for the moment, we proceed to bound the variance of A/" n [c£> 2 ' n ] . In 
what follows, it will be useful to avoid the spectral edges. For this purpose, we let n {x) 
be a smooth, real function such that < 9 n < 1, 

d n (x) = 1, x G [-y/2 + n- 1 / 3 ^ 2 /2, V2 - n" 1 / 3 ^ 2 /2], 

and 

e n {x) = 0, for x G [-V2 + n- 1 ^ 52 ,^ - n - 1/3+52 ] c . 

Then 

Var(AA n [ ¥ p 2 -"']) < 2Var(AU0„< / ? 2 ' n ]) + 2Var(AA n [(l - n ))^ 2 ' n }). 
The rigidity of eigenvalues theorem [H] (Theorem [8] above) , implies that 

with overwhelming probability. Since ||(1 — ^n) < / ?2 ' n ||L cx) < llv 92 '™!!^ 00 ! we have: 

Var(AA n [(l - 6 n ))^ n ]) < \\^ n \\ 2 Lac n 1+452 . 



By the assumption (p G C l l 2+€ and Lemma 3.1, we have, 



\\,„2,n\\2 / ||,„||2 „-l-2e 

It follows that 

provided 2e > 4^. Together with the L°° bound 



II V ,n IU~ < lly 2 ' n | 



L 3 
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this justifies the simplifying assumption, which we make from now on, that p 2 ' n (x) = for 



x 



> ^2 - n - 1 /3+5 2 

Using equations (24), (25), and (26) following the derivation of (23), we write 



(48) 



Var(AT n [^ n ]) = J! \^ n (x)-^ n (y)\ 2 T 2 (x,y)dxdy, 



where 



T 2 = n 2 p n l {x)p n x(y) - n(n - l)p nt2 (x,y), 



and p n ,i, Pn,2 denote the one- and two-point correlation functions of M. For technical 
reasons which will become clear below, we remove a small neighborhood of the diagonal 



x = y in (48). Define 



A = {(x,y) :\x-y\<n 1+5 }. 



We write 



ip 2 ' n (x)-^ n (y)\ 2 T 2 (x,y)dxdy 



_2,n/„.\|2 



<p 2 > n (x)-^ n (y)\ 2 T 2 (x,y)dxdy 



2,n/ \|2 



A c 



h / 1 \^ n {x) - ^ n (y)\ 2 T 2 (x,y)dxdy. 



The "diagonal" term in the previous equation has the following simple estimate: 



(49) 



\<p 2 > n (x)-<f 2 ' n (y)\ 2 T 2 (x,y)dxdy 



n 



\x-y\<n- 1 + s 



f(x, y)(p nj i(x)p n}2 (y) + p 2>n (x, y)) dxdy, 



where \ f(x,y)\ < ||¥> 2 ' n ||i,°o- We cover the diagonal of [—3,3] x [—3,3] by n 1 53 squares 
Sj = Ij x Ij with \Ij\ < n~ 1+<53 , and obtain the bound 



\<p ' \\L°° n 



p n ' 1 (x)p n '\y) dxdy < n- 1 -* ^ (E|A, : A, G Ij 



< n- l - e n l - 5i n 2&3 



< n 



S3— e 



The term involving p 2,n (x,y) in (49) is bounded similarly. 
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Using the determinantal structure, we can express T2 in terms of fC n as 
(50) 

T 2 (x,y) = [ K n {x,x-z)dP( N \z{H)) [ IC n (y, y; z) dP^(z(H)) 



(JC n (%, x) z)K. n (y, y; z) - K n (x, y; z)K n (y, x; z)) dP^ (z{H)) 
JC n (x,x;z)dpW(z(H)) j JC n (y,y;z)dpW(z(H)) 
K n (x, x; z)K n {y, y; z) dP^{z{H)) + / K n (x, y; z)K n (y, x; z) dP^{z{H)). 



We insert the final expression for T2 in (50) into (|48j) and find 
(51) 

Var(AA„[^]) = J! \^ n {x) - ^ n (y)| 2 (| lC n (x,x; z) dP^(z(H)) 

x / K n (y,y;z)dpW(z(H))dxdy- [ K n (x, x; z)JC n (y, y; z) dP^(z(H))) dxdy 



\^ n {x) - V 2 ' n (y)\ 2 J JC n (x, y; z)JC n (y, x; z) dP("> (*(#)) dxdy. 
First, consider the contribution from the final term 

(52) II \^ n (x)-^ n (y)\ 2 j JC n (x,y;z)JC n ( y ,x;z)dpW(z(H))dxdy. 

The saddle point analysis in the next section yields the following bound for the product of 
the kernel evaluated at two off-diagonal points on a set of dP(z)-overwhelming probability. 

(53) \K. n (x, y; z)K. n (y, x;z)\<\x- y\~ 2 
in the region 

\x - y\ > n~ 1+s , \x\,\y\ < V2-n" 1/3+S2 . 



Thus, the second term in ( 50 ) is bounded up to a constant by 

1 



I ||2 „-l-2e 



\x-y\>n~ 1+s \ x V 



2 dxdy < \\^\\ 2 cl/2+e n~ 



for c > 0. We have bounded the contribution to (50) from (51). The error terms resulting 
from sets A of negligible probability are easily controlled by the estimate 



1 2 1 1 2 



\)C n (x, y; z)K n {y, x; z)\ dxdy dP^(z(H)) 



< 11, o 2 11 2 

Wn\\L° 




n 2 (p nA (x)p nA (y) + p 2 ,n(x,y))dxdydP {n) (z(H)). 



Using Lemma 6.4 
(54) 



we can write 

/C n (x, x; z) = n(p sc (x) + ip(x; z)), 
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with |"0(x;2;)| < Cs 1} s 2 n l l 2+&1 for \x\ < \/2 — n V 3 +& and z in a set of overwhelming 
probability. When |x| < \/2 — n _1 / 3+<52 , or |y| < y/2 — n _1 / 3+<52 , we use (54), together with 
the local semicircle law to show that the double integral over x and y of the difference 

J K n {x,x;z)K, n {y,y;z)dP^{z{H))- J JC n (x,x;z)dP^(z(H)) J K n {y,y;z)dP i - n \z(H)) 

is bounded by n 2 n~ 1+252 . In deriving this bound, we must exclude some sets whose com- 
plements have dP( n )(;z(-ff))-overwlielming probability. The corresponding error terms are 
controlled by noting that 



)C n (x, x; z) dx = n, 



whatever the value of z. 

By Fubini's theorem, we obtain a bound of C\\(p 2 ' n \\ 2 ^ 0O n 2 n~ 1+2milx ( S2 ' Sl ^ for first term 
on the right side in (51). Choosing 61,62,63 < e/20, we have established (36). 



7. Saddle point analysis of the determinantal kernel 



In this section, we justify the asymptotic approximations (53) and (54). For this, we use 
2 / — — 

the representation of Kn given by E. Brezin and S. Hikami [4J, and further analyzed by 

K. Johansson [22] : 



(55) 



K l J in (x,y;z) 



e 2n{x 2 -y 2 ) n 



1 



Mf n ,x( w )-fn,y( s )) dwd& 



,,,Jr L w-s 

n 



3=1 



fn,x(w) = 2(w 2 - 2xw) + - V]log(w - zj). 
7 is a rectangular contour 

4 

7 = (J 7i 



i=l 



71 : t -)• —M + it, 
73 = -71 



-ib < t < ib; 



72 : t -> i + ib, 

74 = 72- 



-M <t< M. 



M is chosen such that 7 surrounds all the Z{ 1 < 2 < n, and is the vertical contour 
: t 1— t- it + L, where < M < L is chosen large enough to make Tl and 7 disjoint 
(Figure [J. To prove Lemma |6.4[ we can follow |13| fairly closely. The starting point 
is the following alternate representation of K^^ n , due to Johansson, which removes the 
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-ib 
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74 









Figure 1. Contours 7 and T 



singularity of the integrand: 

(56) K^ 4n (x,y;z) = J j h(s,w)g n (s,w)e n ^ n ^- fn ^ dwds, 



. , 4 / 1^ Zj 

g n {s, w) = -\ w + z-y--y 

z I 4 {w - Zj )(s - Zj ) 



h(s, w 

As noted in [22], we hava^] 



:r- 

e 4n(x-y)w ^ e 4n(x-y)w _ e 4n(x-y)(w-s)^ 

y - x 

fn = fn,x(w). 



g n {s,w) = -f n (w) + 



s s — w 



'Note that there is a typographical error in the corresponding equation in the print version of [23] . 
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The result of the saddle point analysis in [13J is 



1 / T \ e (x 2 - y 2 )/(2a 2 ) 

JC n [x,x+ =- K a J n 



X, X 



n V np sc (x)J n \ np sc (x / 

= p sc {x) S ^{l + o{l)) + 0{n-^). 

TIT 

uniformly for r in a compact set, |x| < y/l — k, and a 2 of order n~ 1+<5 for k > and S > 0. 
We take r = and fix the size a 2 of the Gaussian part to be a 2 = j. On the other hand, we 
want to allow values of x such that \x\ < \/2 — n^ 1 / 3 ^ 2 . This is the statement of Lemma 



6.4 The only significant modification we must make to the argument in |13j lies in the 



localization of the saddle points of f n (s) and f n {w) in (56). In [13], section 3.2, it is shown 
that the equation 



n 

(57) f n {s) = 4( S - x) + - J2 — 







n * — ' s 

3=1 



has two conjugate complex roots s^, which are approximately equal to the roots of 

(58) f'(s) = 4(a - x) + 2(s - ^s 2 - 1). 
The latter are given by 

(59) s ± = -x ± i-\/2-x 2 . 

4 4 

The information we need about s^ is contained in the following 
Lemma 7.1. Let 52 > 0. For any z£Z„ and i£R suc/t t/iai 

\x\ < y/2 - n" 1 / 3 ^ 2 , 

£/ie equation 

_f'n(s) = 

/ias too conjugate roots and s~ = st ■ For large enough n, these roots satisfy: 

(60) | s ±_ s ±|< n -i/2 5 

(61) 9s+ > n" 1 / 6 ^ 2 / 3 . 



Proof. We rewrite (57) as: 



1^1 

s = F n ls) := x — — > . 

w 4n ^ s - z j 

3=1 J 

Following |13j , we show that for large n the mapping F n is a contraction on 

E = {s : |s - s + \ < n~ 1/2 }. 

This will immediately imply the rest of the lemma. It is clear from the explicit expression 
((591) that 9s+ > cn -V6+<5 2 /2 if \ x \ < _ n -i/3+<5 2 . 
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Note that s + is a solution of the fixed point equation 

s = F(s) := x — -(s — \/s 2 — 1) 



The determination of the square root is chosen such that \J s 2 — 1 ~ s at infinity. This 
can be written as y/s 2 — 1 = \/s — lV*? + 1 with the last two instances of %/• denoting 
the principal determination. The condition 9s + > n - 1 / 6 + s 2/3 an( j ^ e \ oca \ semicircle law 
together imply: 

\F n (s)-F(s)\ =0(n- 5 / 6 ) 
for s£S. Moreover, a simple calculation shows that 

\F'(s + )\ < l-Un- 1 ^ 6 ' 2 



for x in the range of interest. By Cauchy's estimate, for s G S, we have 

1 n 1 

i=i v J7 
= F'(s) + 0(n- 2 ^- 6 ^ 3 ) 
1 s 



2 

Thus l-F/^s)! < 1 — 0(n -2 / 3 ) for s G 3. The above implies that 
sup \F n (s) - s + \ = sup \F(s) -s + \+ 0(n" 5/6 ) 

= sup \F(s) - F(s + )\ + 0(n" 5 / 6 ) 

= sup|F'(s)|n- 1 / 2 + 0(n- 5 / 6 ) 

= (l_ cn -l/3+<5/2 )n -l/2 + 0(n -5/6 ) 

< (1 - c'(n))n- 1 / 2 

for some < c'(n) < 1 when is large enough. This establishes the claim. □ 



With (61) in hand, we can follow the procedure in [13], Section 3.4, to localize the 
main contribution to the double integral to a neighborhood of size e = n~ 1+s of the four 



points w = and s = and obtain Lemma 6.4 Some estimates have to be carried out 
differently, because the parameter t appearing in [13J is taken to be small, while it is equal 
to a 2 = \ in our case, but this does not present any serious difficulty. 



We turn to the proof of the off-diagonal estimate ( 53 ) . To obtain a bound on the quantity 

Kn in (x, y)Kn^ n (y, x), we find it convenient to work with the integral representation (55), 
whose form is somewhat more symmetric in x and y. Although the integral is more singular, 
the main contribution to the two integrals is localized away from the singularity for those 
pairs x, y which concern us. Indeed, the next lemma shows that the two saddle points are 
separated by a distance at least a constant times \x — y\: 
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Lemma 7.2. There is a constant C > such that for each n and z G Z n , the solutions 
and s+ of the saddle point equations fn,x( w n) = an d fn,y( s n) = with positive imaginary 
part satisfy: 

\w+-s+\ >C-\x-y\ 
for \x — y\ > n~ 1+s and \x\, \y\ > y/2 — n~ l /^ +&2 . 

Proof. The existence of two con juga te solutions w^, resp. s^ to each of the saddle point 
equations follows as in Lemma 7.1 When \x — y\ > lOn -1 / 2 , the lemma then follows 
immediately from (60). If x and y are closer, we start from the saddle point equations 

1 n 1 

f' n Jw) = A{w - x) + - = 0, 



n * — • w — Zn 



1 " 1 

4(s - y) + - V — 



0. 



n i — • s — Zj 

J= i j 



Subtracting these two equations, we obtain 



1 n 

w-s = x-y + {w-s)-—Y^ 



^ (S - Zj){w - Zj)' 



Rearranging, we find, 
(62) 

d 



7 

(W - S) [ 1- V 7 

I (s - yj )(w 



x-y. 
1 



dt{s-t)(w-t) (s-t) 2 (w-t) (s-t)(w-t) 2 ' 
Thus, by (61), we have 



The assumption z £ Z n allows us to write: 



1 n 

-£ 



1 



n ^ ( s ~ z j)( w ~ Vj) J (s-t)(w-t) 



1 



n(x) 



d 



Pn,l(dt) 
1 



dt (z - t)(w - t) 



dt 



1 



Psc (dt) + ||i/((« - -)(w - -))\\cm ■ 0(n- 1+s ) 



(s - t)(w - t) 

where n(x) = Pi,n(dt). It remains to compute the quantity 



/ 



(s - t)(w - t) 



Psc(dt) 
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For this purpose, it is useful to develop the integrand into partial fractions: 



1 



11 11 

+ 



(s — t)(w — t) w — ss — t s — ww — t 

Computed separately, the two integrals on the right result in the Stieltjes transform of the 
semicircle distribution, evaluated at s, resp. w. Thus we have the identity: 

1 , . , 2 



(s - t){w - t) 



Psc{dt) 



w — s 



- (s — v s 2 — 1 — w + \/ w 2 — 1) 



-1 + 



w + s 



Using (60), we have 



S+)2 - 1 = V( S +)2 - 1 + 



sup 



Vw 2 - 1 + Vs 2 - 1, 

0{n- 1 ' 2 ) 



\s + + a\ 



la^n-^a \y/(s+ + a) 2 - 1| 



By (59): 



|V(s+) 2 -l| = ^[(10x 2 - 18) 2 + 36x 2 (2 - x 2 )] 1/4 . 



It is readily checked that the function x i— >■ (10x 2 — 18) 2 + 36x 2 (2 — x 2 ) is symmetric, 
positive on [0, y/2], and strictly decreasing on that interval. At x = \/2, it takes the value 
4. Thus, we have: 



(63) / - + Sn i = /n/rr/rnri + 

Since w+ = s+ + Oin- 1 / 3 - 6 ^ 2 ), the left side of p3) can be approximated as 



w + + s 



1 



+ 0(n 



-l/3-«5 2 /2 



)• 



The real part of this quantity is bounded above by 3, and is strictly smaller when |x| < \/2 

+ 0(n -l/3-«5 2 / 2) ^ 



(see (77) below). The equation (62) becomes 

y 3 s + 



(w — s) 



2V(s+) 2 - 1 



x- y. 



The result follows by taking absolute values on both sides of the last equation. Moreover, 
if \x\, \y\ < y/2 — £q for a fixed eo > 0, then we also have the reversed inequality: 



\s + — w + \ < C eo \x — y\ 
The final result of this section is the following 



□ 
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Proposition 7.3. Let 8,62 > and z 6 Z n . There is a constant Cs,s 2 such that, for any 
x and y satisfying 

\x — y\ > n~ 1+s , 

\x\,\y\ > y/2-n- 1/3+52 , 

we have the estimate: 

(64) \Ki^(x,y-z)K^(y,x;z)\<C s , S2 \x-y\- 2 . 

Proof. We will show that the product \Kn (x, y; z)\ \Kn 4n (y, x; z)\ is bounded by the 
minimum of two quantities, one of which bounded by the left side of (64) for 

\x-y\> 100n" 1/2 , 

the other being bounded by the left side the equation when 

\x-y\< 100n~ 1/2 . 



In either case, the starting point is the contour integral (55). If, in this integral, we move 

the contour Tl into T^y with < V < M, we obtain: 

(65) 

K^(x,y) = eM f~ V y [ [ -J- e n(/^(«)-/n.,W) dwds + e^-rt shl ^ X ~ »» , 
(v™) 2 J^Jr^w-s ir(x - y) 

The assumption \x — y\ > lOOra -1 / 2 implies, per (60) , |3ft(w + — s^)| > 50n _1//2 . This 
will allow us to choose, for each n, b such as to make the second term in (65) vanish, and 
to extend the integral over 7 to an integral over an unbounded contour 7(6), the union 
of nearly horizontal curves extending from —00 to 00, traversed in opposite directions. 
Indeed, we can find a solution b of 

(66) sin(4n6(x — y)) = 

— 1/2 I I 

within distance ra 100 of any horizontal line (Figure 2h. 

For the purposes of illustration, we assume 3fts+ < 3iw^. We first choose L = $lw+, and 
let the portion of 7 in the upper half-plane be 

71 ■ 1 1 y i^$Sn + t, —00 < t < 3f?s + — e 

ib(n) + t, $twn < t < 00, 



here b(n) is the solution of (66) closest to We choose the part 7* of 71 outside an 

e- neighborhood of s+ and between 3fts + and 5i«;, + to be a smooth curve 7*(i), < t < 1 
with increasing real part, and monotone imaginary part, and with imaginary part always 
lying between 9s+ and b(n). The component 72 of 7 in the lower half plane is obtained by 
reflection about the real axis. We let 

m = |9s + — b(n)\. 

Our assumptions imply that 

k = -»s+| > 50m. 
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w n 



ib 




Figure 2. Contour for the case \x — y\ > lOOn x l 2 



The contribution from neighbourhoods of size e = 1 1 Ks+ — | of the saddle points 
and to the integral 



< /4 >,y) 



e 2n(x-y) 2 n 



(iri) 2 J^j r w-s 



Mfn,x(w)-f n , y (s)) fiyjds 



are now obtained by a standard Laplace approximation, as in [13J. We find that that 
K^ 4n (x, y) is given by the sum of four terms of the form 



a 2n(x 2 -y 2 ) 



2tt 



for all four choices a,b £ {±}. Multiplying this sum by the corresponding approximation 
for Kn^ n (y,x), the real parts of all exponentials involved cancel, since 



(67) 



fn,x( w n) = fn,x( W n) 



fn,y( s n) ~ fn,y{ s n)- 
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A simple calculation, using also the approximations in Lemma 7.1 shows that fn X { w n)i 
fn,y( s n) are bounded below. Finally, by Lemma 7.2 we have 



1 < 1 



K-41 ~ \x-y\ 
for any choice of signs a,b G {±}. It follows that the contribution to 

|^/ 4 ™(x,y)Ky 4 "(y,x)| 

from e-neighborhoods of the saddle points is bounded by \x — y\~ 2 (l+0(n~ 1 / 2 )). It remains 
to estimate the contour integral Kn^ 4n (x,y) for s, w at distance greater than e from w^, 
• We will show that the integrals over these portions of the contours give a contribution 
bounded up to a constant factor by 

1 C nm n .x(w£)-Wn.v(sj)) _ 



\x - y\ 



Using the symmetry (67), the real parts of the exponential factor will cancel with the 



corresponding factors from Kn n (y,x). Define the set 

n = {s = a + i/3 :\a- Ks+I > e > \P ~ % s t\ < 50 -1 e}. 

The part of the contour 71 outside of a ball of radius e around s+ lies entirely in Q. For 
s G fl, we have: 

M[fn,y(a + ip) - f n , y (s+)} > c(q - K S +) 2 

for some constant c > 0. A computation using f'(s + ) = shows that for a + i(3 6 0, 

mf(a + ip) > c(a-Ks+). 

Thus 

d x Uf n , y {a + > Sft/'(a + i/3) - 0(n- 5 / 6 ) 

> »/'(«+) - 0(n- 2 / 3 ) 

> c(a - 

>c'(a-Ks+), a + i|3ea 

Integrating then gives the desired bound. We also get a similar lower bound along the 
contour V (again, see [E]). When we integrate outside the saddle point, one of the terms 
we have to deal with is of the form 

dw [ dse 11 ^*^- 1 "^— 1 —, 



r Jj w - s 

where 7 denotes the part of the contour 71 outside a neighborhood of of size e. We will 
derive the bound for the above integral, the case where the w lies outside an e-neighborhood 
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of the saddle point is essentially the same. Using the bound derived above, we have 

1 



dw / dse n ( fn '*W- fn ^ - 



w — s 



1 



w — s 

Split 7 into % and where j Q is the part outside of the (integrable) singularity, and % 
is the part within j^q\^(w^ — 4)1 of the singularity. The above integral can be expressed 
as: 

(69) 



w — s 



1 

u; — s 



For the first term (68), we use that , j_ . < "V 00 + on the contour, and 

I I \ji\'W n s n )\ 

,-cn(5R S -5Rs+) 2 ds < n ~l/2_ 



To 



We obtain an additional factor of n -1 / 2 for T from either the e- neighborhood saddle point 
io„ or the gaussian decay of the w integral away from • We conclude that the term ( 68 ) 
is bounded by 

lOOn" 1 ^^ - s +)|- 1 e" ffi (/n.-(^)-/»,«( s «)) < l00n _1 |x - y^e 7 ^^'^-^^. 

For the second term (69), notice that since |3?s — < j^q\^{w^ — 4)1) and idGT has 
constant real part, we have 

|5fts - sRs+| = |Ru;+ - 3fo+ + Rs - &w\ 



>i|»(«,+ _4 



l 

Too 



!»(«;+ -4)1 



Therefore, we have |5fts — ;R4| 2 > j|5ft(/i4 — 4)| 2 - The singularity behaves like 1/r 

in two dimensions, therefore the second term is bounded by 



From the fact that 



x 2 e -nx < 



for x G M, we conclude the estimate: 

e 



9 f to(/n.s(Wn)-/n.p(«n)) ^ j^-f _ g + | e C7l|5R(«)+ -s+ ) | 2 < e r*(/„, a; (w^)-/n,„(sJ)) n -l |S|(J( U ;+ — 4)| _1 



< e nS(/„,,(«i{)-/n 1 ii(»»)) n -lL _ y|-l > 
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In the case \x — y\ < 100n -1 / 2 , we do not attempt to set the sine term in (65) to zero. 
We choose 7 to consist of two horizontal lines through s+ and s~, respectively, traversed 
with opposite orientations (Figure [3]). We can then reproduce the saddle point analysis 



-b 



Figure 3. Contour for the case \x — y\ < lOOn 1 I 2 

for Kn An (x,y) and Kn in {y : x). When we multiply the resulting approximations, we find 
16 terms similar to the ones encountered above, each of which is bounded, up to constant 

factors by , — . There are also 8 terms of the form 

\x-yr 

Cn ,b Sin(2n5s+(g - y)) g 4n(y-a!)«u.+ e Ti(/n, v (^)-/n,«i«)) 

<-4 n(x-y) 

where a, b E {±} and c^ b are bounded independently of n. All these terms are bounded 
by \x-y\~ 1 U P ^° a cons tant factor. To see this, it clearly suffices to show that 

(70) n • (4(y - s)^wt + ®-fn,y{s h n ) - R/^K)) 

is bounded by a constant independent of n, for any choice of signs a, b £ {±}- Notice the 
identities 

fn,x(w) = fn,x(w), 
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and 

so that we have 



fn,y{s) — fn,y{s), 



^fnA w n) = ^fn,x( w n) 

Thus it suffices to deal with the case a = b = + in (ffol) . Write 



1 n ( — \ 
fn, y (s) - f niX (w) = 2(s 2 - 2ys - w 2 + 2xw) + - V log f 1 + 

J=l v J' 

1 n 1 

= 2(s 2 - 2ys - w 2 + 2xu>) + (s - w) ■ - ^ 



n » — * w — Zj 

3 =l J 



(s — w) 
2n 



- E 7-^2 + 1 E E(- 1 ) fc ~ 1 ( — 1 

j=l v ■ ?/ j=l k>3 



When w = w£, we can use the saddle point equation (57) to replace ^ J2]=i w - y - m the 
above equation by 

4x — 4u>+ . 

Thus: 
(71) 

fn,y( s n) ~ fn,x{ W n) = 2 (4 ~ *4 f + 4 ~ 2/)4 



(4 - W n ) : 

2n 



" 1 1 n /«+ 71,+ V 



3=L j, j. 

Taking the real part of the first line on the right, we find 

(72) 2 ((»(,+) - K«)) 2 - (9(4) - ^(wt)) 2 ) + 4(x - y)R.+ 

Note that the term involving imaginary parts is always negative. We will find below that 
it cancels another, positive and real, term in the sum. Using (60), we have: 

»(*+) - M(w+) = l( x - y ) + 0(n-V3) = O^- 1 / 2 ), 

9(4) - 9f(w+) = -(y/2-y 2 - V2 ~ x 2 ) + Oin- 1 / 2 ) = O^- 1 / 3 ^ 2 / 2 ). 



By (72), we see that we can estimate (70) provided that we can control the real parts of 
all terms appearing on the second line of (71). Since we have 

| s + _ ^+1 < „-l/3-fa/2 



W 



■+- Zj \ >%w+ yn- 1 ^/ 2 , 
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we can easily estimate the sum 



n j=i k>3 \w n -z 3 J 



< 



-l/3-<S 2 \Jfe 



k>3 



< n -l-35 2 _ 



_ n 

It remains only to deal with the term 

(4 ~ wj? y 1 

2n j^ {w +_ z . ) 2- 

The potentially dangerous terms comprise the real part of the product above: 



(73) 
(74) 



1 1 n 1 

- ((»(*+) - M(w+)f - (9(4) - 9(4)) 2 ) »- E — r 

i n i 

n 7=1 l^n ~~ z jJ 



As in the proof of Lemma 7.1, we have the approximation 



1 1 



(75) 



-1 + 



-1 + 



-71 



:^) 2 - 1 

v^) 2 - 1 



+ 0{n-^- s ^) 
+ 0(n-V2), 



where s + is as in ( 59 ) . We examine the imaginary and real parts of the quantity ( 75 ) . By 

a/(s+) 2 - 1 = 3s+ - 2x. 



(58), we have: 



(76) 
Thus, we have: 

(77) 



3a; 2 + 3V2 - x 2 9xy/2 - x 2 
i- 



v /( s +)2 _ i 3s+-2x x 2 + 9{2-x 2 ) x 2 + 9(2-x 2 ) 

The real part is a symmetric, positive function which increases from to 3 on the interval 
[0, s/2]. The imaginary part is an anti-symmetric, negative for < x < \[2 and zero at the 
endpoints. 



For (73), we find: 

(9(4) - Q(*4)) 2 



3x 2 + 3V2^ \ _! 
x 2 + 9(2-x 2 ) - 1 + ° (n } 
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The factor in parentheses is at most 2, and so the sum of the first term in the previous 
equation with the term — 2(Q(s + ) — Q(w^) 2 from (72) is at most zero. For the term (74), 
we find: 

(K4 - »u;+)(9te+ - 9to+)9f- £ ; + 1 , 2 = -j(K*t ~ ^+)(^^ - >/2=tf) 

n ^ (w£ - Zj) 1 I 



9xV2 - x 2 



+ 0(0 



x 2 + 9(2 -x 2 ) 

v / 2^ 2 1 



-$ls+-ttw+){y-x){y + x) 



+ 0{n- 1 ) 



For 2 — y 2 = o(\x — y\), the ratio of radicals is bounded, and so the last quantity is 0{n~ l ). 



All the terms in (70) are now accounted for, and the argument is closed. □ 

8. CLT FOR C l ~ e : WlGNER MATRICES 

In this section, we prove the CLT for general Wigner matrices under matching moments 
conditions at the C 1_e level of regularity. The proof follows similarly to that of the previous 
section. Writing 

\fc<(l+5)logn fc>(l+<5)logn/ 
= <p 1 ' n + <f?' n , 



where 5 is as in Proposition 6.1 The variance of the first term is bounded as in Section 
[6j we now compare M with a GUE matrix with entries of complex variance 1 rather than 
l/y/2 times a GUE matrix. For the second term, we use the trivial bound: 

Var(AA n [^])<n 2 ||^|| Loo 

<n 2 £ 2- 2 ^-)|l^l^ 

k>(l+8) logn 

Picking e sufficiently small, the above quantity is o(l). As in previous cases, this bound 
for the variance implies the CLT. 
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